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If} , We study the resonant divergences that occur in quantum scattering cross- 

sections in the presence of a strong external magnetic field. We demonstrate 
that all such divergences may be eliminated by introducing radiative correc- 
tions to the leading-order scattering amplitudes. These corrections impose a 
choice of basis states that must be used in scattering calculations: electron 



5/5 , states must diagonalize the mass operator, while photon states must diago- 

nalize the polarization operator. The radiative corrections introduce natural 
line- widths into the energy denominators of all propagators, as well as into 
^ ■ the time-development exponentials of all scattering states corresponding to 

external lines. Since initial and final scattering states may now decay, it is 
logically necessary to compute scattering amplitudes for a finite time-lapse 
between the preparation of the initial state and the measurement of the final 
state. Strict energy conservation, which appeared in previous formulations 
of the theory, must thus be abandoned. We exhibit the generic formulae for 
the scattering cross-sections in two useful limits, corresponding to the cases 
where either the initial states or the final states are stable, and discuss the 
application of the general formula when neither of these limits applies. 

12.20.-m,11.10.st,97.60.Jd 

Typeset using REVTgX 



*E-Mail Address: carlo@twinkie.gsfc.nasa.gov 



1 



I. INTRODUCTION 



Astrophysicists have had a long-standing interest in the physics of elementary processes in 
super-strong magnetic fields, with field strengths B > 10 12 G. The cyclotron lines observed in 
the spectra of Her X-l |l| and of 4U 0115+63 0] as well as in many other X-ray pulsars have 
energy centers which correspond to field intensities in this range. There is also evidence for 
such field strengths in the spin-down rates of radio pulsars. If the spin-down is attributed to 
energy loss to electromagnetic radiation from a spinning magnetic dipole, many observations 
are consistent with field strengths of the order of 10 12 -10 13 G, with some pulsar field strengths 
well in excess of even 10 13 G In addition, there is tantalizing evidence for cyclotron 

lines in the spectra of gamma-ray bursts seen by the gamma-ray burst detector aboard the 
GINGA satellite ||, with line center energies consistent with field intensities of order 10 12 G. 

The association of Soft Gamma Repeaters with supernova remnants provides indirect 
evidence for even stronger fields. If the 8 s periodicity of the March 5 1979 event is identified 
with the rotation period of the neutron star, the known age of the N49 remnant may be 
used to estimate that the field strength is approximately 6 x 10 14 G @. Moreover, such a 
strong field could help resolve the puzzle of how the March 5 1979 event could ostensibly 
have been so extravagantly in excess of the Eddington limit (L « 10 4 LEdd)j by suppressing 
the Thomson cross-section for photons propagating nearly parallel to the field lines J7J. 

In fields such as these, comparable in strength to the critical field strength B c = 
m 2 c 3 /eh = 4.414 x 10 13 G, all calculations of elementary processes must be carried out 
using Quantum Electrodynamics. There have been many such calculations over the past 
two decades, covering topics such as cyclotron absorption ||, cyclotron decay [PJIJ, sin- 
gle photon pair production [11], 12], pair annihilation to a single photon [13—16 1 , Compton 



scattering |T7|-pO|, two photon pair production two photon pair annihilation [|l"3"l , |l4| 
e~e~ scattering |EJ, and several more. All these processes have very different behavior from 
their B = counterparts (if those counterparts are even possible), on account of the peculiar 
kinematics, as well as the discrete electronic states (Landau levels) associated with a uniform 
external magnetic field. 

These calculations have always been carried out in the Furry picture, in close analogy to 
the B = Feynman rules. The free space electron propagator is replaced by a propagator 
which is a Green's function for the Dirac equation in the external field, and the external 
fermion lines are represented by solutions of that equation. The results have often been 
interesting and useful, but they have not been uniformly satisfactory. The leading order 
calculation of resonant Compton scattering yields results which are divergent at the cyclotron 
resonances , evidently because to this order the theory makes no provision for natural 

line width. The line width may be included "by hand" in the results, making the cyclotron 
resonances finite [19,2(j]. However, there are other such "resonant divergences" in Compton 
scattering, which have nothing to do with the cyclotron resonances and which are not nearly 
as tractable |17| , |l8| . In fact, one such resonance, which occurs exactly at the threshold 
where the initial photon may pair-produce, is responsible for making the total Compton 
cross-section divergent everywhere above this threshold. 

In fact, the theory of elementary processes in external magnetic fields is plagued with such 
divergences. With a little practice, it is not hard to discover divergences resembling these 
in every single process of second or higher order. Clearly, this is a troublesome development 
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that casts a shadow upon the entire undertaking. 

These resonant divergences occur because the kinematics of these processes allow inter- 
mediate "virtual states" to be real — that is, on-shell. Often, the on-shell intermediate state 
is an excited state, such as an electron in an excited Landau level, or a photon above the 
one photon pair-production threshold. In this case there is an associated decay width that 
may be pressed into service to control the divergence. The Weisskopf-Wigner broadening 
prescription, 

E^E-^iT, (1.1) 

when applied to the energy denominators in the propagator, pushes the poles in the propaga- 
tor off the real axis, so that while the intermediate states may still be on-shell the propagator 
no longer diverges there. This is in fact the approach that has been adopted for Compton 
scattering [19,20|23|, and which ascribes the natural line width to the cyclotron resonances. 



However, there are circumstances in which a stable on-shell intermediate state may be 
produced. Such states are not attended by a decay width, so the associated divergence may 
not be reined in as before. 

It should be pointed out that these divergences are entirely unrelated to the notorious 
ultraviolet divergences of QED. They are not the consequence of improper manipulation of 
field-theoretic distributions; rather, they occur whenever the circumstances of the elementary 
process permit a kinematically accessible on-shell intermediate state (KAOSIS). It is easy 
to recognize when a KAOSIS is permitted. For example, a second-order process will allow 
one if it may be viewed as a succession of two real first-order processes. Thus the KAOSIS 
corresponding to cyclotron resonance occurs because the process may be viewed as a real 
cyclotron absorption followed by a real cyclotron emission. Similarly, the second "disastrous" 
resonance in Compton scattering is due to a KAOSIS that corresponds to the initial photon 
undergoing a real decay to a pair, followed by the resulting positron annihilating with the 
initial electron to produce the final photon. The reason this KAOSIS is catastrophic is that 
the intermediate positron may be in the Landau ground state, so that no decay width is 
available to restrain the divergence. 

At the same time, there is a second defect of the theory which so far has not received 
recognition as a problem. The calculation of S-matrix elements as outlined above always 
results in a ^-function which enforces strict energy conservation between initial and final 
states. This remains true even if some of these scattering states are unstable. But this is 
not physically sensible; the energy of an unstable state is only known to within its decay 
width, so that it is ludicrous to demand strict energy conservation for such transitions. 
Nevertheless, the current theory does so irrespective of whether or not the states are stable. 
As an example, the calculation of the decay of an electron in an excited Landau state 
produces the result that the emitted cyclotron photon is monochromatic, rather than having 
the Lorentzian line shape characteristic of resonant decay P,[10 . 



The difficulties of resonant divergence and of spurious energy conservation are related. 
Briefly, a stable KAOSIS may only occur if some of the particles in the initial and final states 
are themselves unstable, and it is their decay widths that restrain the divergence. However, 
the introduction of these decay widths smears out their energy, so that energy conservation 
(which was a consequence of their assumed eternal duration) no longer obtains. 
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Thus it appears that we must modify the current theory somehow if we are to circumvent 
these unphysical features. That modification is the purpose of this work. We demonstrate 
below that radiative corrections modify the propagators and the scattering states by in- 
troducing their respective decay widths into the S-matrix elements. In this sense, we are 
extending the results of Graziani [p]], who carried out this program for the electron propaga- 
tor only. The corrected "dressed" states and propagators result in scattering cross-sections 
that are always finite. 

In Sec. |T| we discuss the role of the bare propagators in producing resonant divergences. 
In Sec. [TIT] we review the work of reference |23| on the electron propagator, which we extend 
to the photon propagator in Sec. |V[ We exhibit the corrections to the scattering states in 
Sec. [V[ In Sec. [VI] we derive the modification of the S-matrix elements, and exhibit two 
useful limits: the absorption limit and the emission limit, wherein the initial and final states 
are stable, respectively. Finally, in Sec. [VI I| we discuss the applicability of these results and 
the extent of the modification from the "standard" theory. 



II. RESONANT DIVERGENCES AND PROPAGATORS 

We begin with a discussion of the role that the bare electron and photon propagators 
have in producing the resonant divergences. In this section and throughout the rest of this 
work we use the metric signature [+1,-1,— 1,-1]. 

The bare electron propagator is simply expressed in terms of the fermion one-particle 
states. We work in the Furry picture, so these states are represented by solutions w2 of 
the Landau-Dirac (L-D) equation, that is, the Dirac equation in the presence of a classical, 
uniform, external magnetic field: 

[r^-eA^-m]^ = 0. (2.1) 

Here, i] = ±1 refers to whether the solution has positive or negative energy, while A is 
the set of quantum numbers that specify the state — Landau level number, x 3 -momentum, 
spin, and an orbit center coordinate. We choose the gauge A = 0, A = Bx 1 e2 for the 
external magnetic field, which is directed in the e3 direction. We also choose the orbit center 
coordinate a to be the x 1 Cartesian coordinate of the orbit center. With these choices, the 
functional dependence of the solutions of Eq. ( [2.1| ) is that of a plane wave in the x 2 and 
x 3 coordinates, while the x 1 dependence is that of a one-dimensional harmonic oscillator 
eigenfunction centered about x 1 = a p5|j2"6| . 



It is convenient to separate out the time dependence of ^a- 

^ ) (x) = e*°#(x), (2.2) 

where 0^(x) is an eigenstate of the Landau-Dirac Hamiltonian with quantum numbers A 
and eigenvalue Ea- Ea is given by Ea = [m 2 + (p 3 ) 2 + 2mNuoB\ 1 ^ 2 ■, where N is the Landau 
level number, p 3 is the momentum in the x 3 direction, and ujb is the Larmor frequency. The 
bare electron propagator Gb(x,v) may be represented in terms of the (jy^\ 



^ J 2tt p° — T)(Ea — i0) 
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It is easy to verify that the expression in Eq. (|2.3|) is a Green's functions for Eq. (|2.1|) with 



the required boundary conditions fl24|| . We write Gs(x,y) rather than Gb(x — y) because 
neither Eq. Q2.1| ) nor its Green's functions are translationally invariant — they are invariant 
under combinations of translations and gauge transformations. 

The bare electron propagator enters the leading order calculation of such processes as 
Compton scattering, pair annihilation to two photons, and two photon pair production. 
The S-matrix elements for these processes are obtained by sandwiching Gs{x,y) between 
appropriate fermion and photon states and integrating over x and y. These processes all 
exhibit resonant divergences. The proximate cause of those divergences is evidently the 
energy denominator in Eq. ( |2.3| ). The convolution of the propagator with the initial and 
final states fixes the values of p° and p 3 at the total energy and x 3 -momentum of the states, 
respectively. Thus, the residual energy-related degree of freedom in the summation over 
intermediate states is the Landau level number N, a non-negative integer. It follows that 
the energy denominator in Eq. ( |2.3| ) may be zero if the energy and x 3 -momentum of the 
scattering states are suitably tuned to one of the Landau levels. If this happens a resonant 
divergence occurs. 

The bare photon propagator has the usual "transverse" expression 

/jAh, -ik(x-y) 
^(-^ + ^ rA2) ___. (2 .4) 

It enters the calculation of S-matrix elements for processes such as electron-electron scat- 



tering and electron-positron scattering. These processes have total cross-sections which 
exhibit resonant divergences. This may seem surprising at first, since the S-matrix elements 
for these processes are entirely finite, even when a KAOSIS is present. There is a key differ- 
ence between the photon and electron propagators: In the electron case the on-shell energy 
contains a dependence on the Landau level quantum number, a discrete degree of freedom. 
Thus the resonances are "spaced-out" and well-separated. On the other hand, the photon 
on-shell energy is entirely dependent on the purely continuous degrees of freedom k, and 
the "sum over intermediate states" is really an integral. While this integral will undoubt- 
edly encounter any KAOSIS singularity, the small imaginary part in the denominator of 
the propagator provides a perfectly straightforward prescription for circumventing the pole. 
Thus the S-matrix elements are finite for these processes. Nevertheless, whenever an on-shell 
intermediate state is kinematically accessible, the total cross-sections diverge. What is going 

The divergence is evident in the expressions in Appendix C of Langer E2] for the e-e 



on? 



scattering cross-section. Langer performed the integrals over the intermediate states simul- 
taneously with those over final states and the average over initial states, in order to take 
advantage of several simplifications that ensue. Unfortunately, these manipulations obscure 
the source of the divergences. In order to shed some light on the situation, we have calcu- 
lated cross-sections for e-e and e + -e~ scattering by computing the S-matrix elements before 
summing over final and averaging over initial states. 

What we have found can be best illustrated by considering the specific example of e + - 
e~ scattering. As depicted in Fig. [I], there are two relevant Feynman diagrams. We have 
confirmed that in each case the S-matrix elements are finite even when there is a KAOSIS. 
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When a KAOSIS exists, however, there is a divergence resulting from the sum over final 
states. In the diagram of Fig. 0(a), the divergence arises from the integral over the variable 
s = — aj , where and aj are the orbit center x 1 -coordinates of the final positron and 
electron, respectively. In the diagram of Fig. 0(b), the divergence is due to the integral over 
the variable t = aj — a~l , where a~ is the orbit center ar-coordinate of the initial electron. 

We see, then, that the divergence manifests itself as an infinite range of interaction. 
When a KAOSIS is present, the intermediate photon may travel an arbitrarily large distance 
between its emission and its absorption without a reduction in amplitude. Consequently 
there is no spatial cutoff in the interaction, and the resulting integral over orbit center 
separation diverges. 

The situation may be further clarified by explicitly computing the interaction range in 
the x l direction. When we calculate an S-matrix element for these processes, we integrate 
the bare photon propagator Db{x — y)^ v multiplied by two fermion currents, one for each 
vertex. Since our states are chosen to be plane waves in the x 2 -x 3 direction, the currents are 
also plane waves in x 2 and x 3 , and the integrals over x°, x 2 , and x 3 are tantamount to Fourier 
transforms of the photon propagator in those coordinates. The remaining integral over x l 
folds together two (in general well-separated) current functions ^(x 1 )^, j y {y l ) u with the 
transformed propagator fs^x 1 — y 1 ; k°, k 2 , k 3 )g^ u . We may easily compute fsix 1 ; k°, k 2 , k 3 ). 
If we define p = (k ) 2 — (k 2 ) 2 — (k 3 ) 2 = ±q 2 , with q > 0, we find 

{_J_ e iq\x 1 \ :f > Q 
Ll. < 2 ' 5 > 

Thus, if p > (the condition for the existence of a KAOSIS), the effective interaction range 
is infinite, whereas it is of finite range ~ \ jq otherwise. 

The situation is reminiscent of Coulomb scattering. If a wave packet scatters with large 
impact parameter from the center of the Coulomb potential, the momentum transfer is very 
low, and thus the exchanged photon is very close to its mass-shell (that is, its k vector is 
very near the apex of the light cone). This is the reason that the total Coulomb cross-section 
is divergent — the extreme infrared photons, which are nearly on-shell, give the interaction 
an infinite range. There is one noteworthy difference here, however: for the processes we 
consider, the virtual photons do not need to be infinitely soft to be on-shell, so that these 
resonant divergences are in no sense infra-red divergences. 

As a consequence of this analysis, the cure for photon-propagator-related resonant di- 
vergences will necessarily have a slightly different mathematical character from that for the 
divergences that arise in connection with the electron propagator. While in the latter case 
we need to show how the previously diverging amplitudes may be made to converge, in the 
former case we must show how the previously long-range "effective interaction" may have 
its range curtailed. 



III. THE DRESSED ELECTRON PROPAGATOR 

The necessity of complexifying the energy denominator of the electron propagator so 
as to obtain the Weisskopf-Wigner line shape has been recognized for some time []H||2"0 
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although these authors lacked a formal justification for their ad hoc complexification of the 
energy, given in Eq. (|1.1|) . Graziani |23 provided this justification by considering radiative 



corrections to the electron propagator. As it turns out, the same procedure may be applied to 
the photon propagator as well. Quite generally, one finds that the operator which represents 
the self-interaction (the mass operator E in the case of the electron, the polarization operator 
IT in the case of the photon) singles out those solutions of the wave equation in which it 
is diagonal. The dressed propagator may then be expressed as a sum over all such states, 
with the energy in the denominator acquiring an imaginary part given by Eq. (|1 . 1| ) . We 
review the discussion of the electron propagator from Ref. pfl , and extend it to the photon 
propagator in the next section. 

The dressed propagator G(x, y) is related to the bare propagator and the self-energy 
operator £(x, y) by the Dyson equation G = Gb + Gb ■ E • G, where the dot denotes four- 
dimensional convolution as well as spinorial matrix multiplication. It follows that G(x,y) is 
a Green's function for the dressed L-D operator 

[ 7 "(id„ - eA^) - m - £] • = 0. (3.1) 

Now, let 9$ A (x) = e~ ip ° x °^\x.). It may be shown [0,0,0 that 6 diagonalizes E 
if the <fr A are simultaneous eigenstates of the Landau-Dirac Hamiltonian and of the x 3 - 



component of the magnetic moment operator of Sokolov and Ternov f28fl . This condition 
defines what is meant by "spin up" and "spin down", in the absence of the tools provided 
by Poincare group representation theory (the physical system no longer has full Poincare 
invariance). If this condition is satisfied, and only in this case, G may be represented in a 
form analogous to Eq. (|2.3|): 



y > ~ J 2vr f - V (E A - iO) - E(pP, A, rt) ' ^ 

A,r) 

where E(p° , A,rj) is the diagonal matrix element of E in the state &pp A - It is in general a 
complex number, in contradistinction to the case where the external field strength is zero, 
where the diagonal matrix elements of the self-energy operator are real. The real part of 
Ti(p°,A, rj) merely yields a small shift in the energy of the state, while the imaginary part 
provides a line-width to the otherwise divergent resonant energy denominator. This justifies 
neglecting the real shift while preserving the imaginary width. It was shown explicitly in 



Ref. [23[ that in this approximation, the pole in the propagator corresponding to state A, 



rj is located at p° = r](E A — ^T A ), where is just the Weisskopf-Wigner decay rate of the 
state A, computed to leading order: 

r\4 = e 2 J ^^J2J2\ T ^e)\ 2 (2n)5(E A -E B ~uj k ). (3.3) 

Here, T A b{Il, e) is the interaction matrix element for the transition from the state A to the 
state B with the emission of a photon with wave vector k and polarization e. It follows that 
the dressed propagator may be expressed as 
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A,T) J 



-ip (x -y°)_ 



(.V), 



(3.4) 
(3.5) 



to first order in e 2 . 

It appears from these equations that the prescription of Eq. ( |1 . 1|) is in fact rigorously 
justifiable. We wish to emphasize, however, that Eqs. ( |3.4j )-( |3~5D are only correct when the 
states <f>jp are chosen so that the diagonalize the self-energy operator, or equivalently, 

so that the (j) A diagonalize the x 3 -component of the Sokolov-Ternov magnetic moment 
operator. The states of Sokolov and Ternov 



2c , and those of Herold, Ruder, and Wunner 



T0|1 satisfy these conditions, and their use leads to correct expressions for the scattering 



cross-sections. The states of Johnson and Lippmann [25], which have gained some currency 
in the literature, do not satisfy the required conditions p6| . In particular, as discussed in 
23] , the use of Johnson-Lippmann states in the computation of cyclotron scattering cross- 



sections can lead to relative errors of order 45% at the first cyclotron harmonic, depending 
on the field strength. 



IV. THE DRESSED PHOTON PROPAGATOR 

The procedure for the photon propagator is entirely analogous to the one followed for 
the electron propagator. The self-interaction of the Maxwell field is represented by the 
polarization operator U(x — y)^ u , which satisfies the condition of gauge invariance, -^ll(x — 
y)^ v = 0. The dressed photon propagator D(x — y) pu is obtained from the bare propagator 
Db{x — y)^ and EE (a; — y)^ by solving the Dyson equation D = Db + Db'^- D. In order to 
accomplish this, it is necessary to find the polarization states which diagonalize IP 1 ', which 
are analogous to the spinor states that diagonalize S. These polarization states were found 
by Batalin and Shabad f29|, and written explicitly for the case of a uniform magnetic field 



by Shabad J30|. For a photon with 4- momentum k propagating in a uniform magnetic field 



B = Be 3 we have the following three (unnormalized) polarization vectors: 

b1 = k 2 ( ei r-k\e 2 r, (4.1a) 



6jf = fc 3 (e r + fc°(e 3 r, (4.1b) 

b1 = (k v k v )k^-(k ± v k ±1/ )k", (4.1c) 

where k±^ = fc 1 (ei) M + fc 2 (e2) At , and (e p ) M is a unit vector in the x p direction. These modes 
diagonalize the Fourier-transformed polarization operator U(k) fJ-l/ = J d A xU(x) tJ-u e lkx . There 
are three of them, because by gauge invariance n satisfies k f JJ(k)^ u = 0, so that the fourth 
mode is just k, and the eigenvalue of n that corresponds to it is zero. The mode b\\ is 
so labeled because on shell it differs from the usual "parallel" polarization 3-vector by an 
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inessential multiple of k, while b± is just the usual "perpendicular" mode. The vector 
represents a longitudinal mode, which on shell is proportional to k. The b in Eqs. (|4.1|) are 
orthogonal to each other, and to k. 

Using these modes, the transverse photon propagator may be expressed as follows: 

D(x _ y v = - [ ^ e -iHx-v)X^h3. I (4 2) 

1 V) J (2vr)4 e ^b r bj (A;T-^ 2 + n(P,k,i)' K } 

where 

n(*°,k,j) = n(*)^6 i 7(6 i -bj). (4.3) 

To first order in e 2 , the pole in Eq. ( |4.2|) is located at (A; ) 2 = cj fc 2 — n(co> fc ,k, j), so that to 
this order Eq. (|4.2|) may be written 



V[x ~ y) = ~JWf e (fcO)2- Wfc 2 + nK,k,j)- ( j 

Note that when the light-cone condition = ujk is satisfied, both 6j_ and b\\ are space-like 
(6 ■ b < 0) while the longitudinal mode bi is light-like. It follows that H(u>k, k, L) = 0, so we 
only need compute H(u)k, k, j) for j =_L, ||. 

The (unrenormalized) leading-order expression for II(x — ?/)'"' is 

IL(x - ^ = -ze 2 Tr(G B (2/, 1)7^(1, y)-/*). (4.5) 

Following the analogy to the case of the electron propagator, we calculate the imaginary part 
of the diagonal matrix elements of IP 1 ', while neglecting the real part. For this purpose, Eq. 
( |4.5|) is entirely adequate, even though it is not renormalized. The renormalization counter- 
terms which are to be subtracted from the diagonal matrix elements of IP^ are purely real, 
so that the imaginary parts are unaffected by renormalization. 

This expression for IP 1 ' is in fact translationally invariant, even though Gb is not. The 
translational invariance may be established using the "translation+gauge transformation" 
invariance of Gb alluded to in the previous subsection [pi]] . 

Substituting Eq. (p.3[) into Eq. (|4.5|), after some manipulation we obtain 



uu - «)»» - -e 2 V (*t e -w-V) [tfW^Wll^'iyl^Ry)] , 4 6) 

[ V) ~ A % J 2vr rrf> + E A + E A ,-iO ■ (U)) 



We now Fourier transform this equation. Taking the implicit translation invariance into 
account, we obtain 

Il(kr^u = L-'T- 1 J d A xd A ye lk ^ y) Ii{x - yT ^ v 



I HrW I 2 
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where 



x e 



ikx 



(4.8) 



When e M is a normalized polarization vector, J^(k) M e M is just the interaction matrix element 
for a transition from a photon with wave vector k and polarization e to a pair with quantum 
numbers A' A. The imaginary part of Eq. ( [4.7[ ) is 

lm[U(kre,e u ] = -2u k J2 Uaa'^Y^ ' ^ (&A + E A > ~ \k°\) , (4.9) 



A.A' 



where we have used the identity J2 AA , |J^,(k) M e M = J2aa> ^aa'O^Y 6 
of the parity invariance of Eq. ( |2.1| ) . 

Comparing with Eq. ( |4.3j ), we see that to obtain the imaginary part of H(w k , k, j), we 

may impose the light-cone condition k° = u k and let e M = (e,-) M = (bj)^/\bj ■ bj] 1 / 2 



r(-) 



, a consequence 



in 



Eq. ( [4.9|) , keeping in mind that bj ■ bj < 0. The result is 

Im[II(u; fc , k, j)} = 2u k J2 Uaa^Y^)^ ' ^ {E A + £ A , - u k ) 



--ui k 



A.A' 



= 2^xr(k,j)/2. (4.10) 

Clearly, T(k, j) is just the Weisskopf-Wigner decay rate of the photon state (k, j). The energy 
denominator in Eq. is thus (k ) 2 -u k 2 + 2iuj k T(k,j)/2 « (A; ) 2 - - «r(k, j)/ 2 ) 2 , to 
first order in e 2 . Consequently, the dressed photon propagator is 



D(x-y) 



d 4 k 
WY 
d 3 k 
(2n) 3 2u k 

Y.Y. ;-:,<■.» 

j r?=±l 3 3 



bj ■ bj (k' 



0\2 



(Lu k -iT(k,j)/2) 2 



(4.11) 



x ° -y )] e -^K-«r(k,j)/2)(x°-y ) e ir,k-(x-y)_ 



(4.12) 



We see that the prescription of Eq. ( |1 . 1| ) continues to hold in the case of the photon field. 
It should be emphasized, however, that it is essential that the polarization modes given in 
Eqs. (|4.1|) be used in the expression for the propagator in order for that expression to be 
correct. 

Recalling the discussion at the end of Sec. U we investigate the range of the "dressed" 
interaction by computing the partial Fourier transform of the propagator in Eq. ( |4.11| ) with 
respect to x°, x 2 , and a; 3 . Assuming the presence of a KAOSIS (so that q 2 = (k ) 2 — (k 2 ) 2 — 
(k 3 ) 2 > 0), we find for the transformed propagator f(x 1 ; k°, k 2 , k 3 ), 



f{x ; k°, k , k 3 ) 



2Z(q) 



Z(q) = -j= [yju* + q 2 + 



i\f u 2 — q 2 



u 



g4 + (fc0)2r(k,j) 



(4.13) 



10 



where T(k,j) is evaluated at the point k 1 = q. Since the imaginary part of Z(q) is positive, 
we find that the presence of a non-zero decay rate r(k, j) has cut off the interaction range. 

Note that Z(q) — > q as T — > 0, so that if the KAOSIS is below the one photon pair- 
production threshold the interaction range is not cut off (in fact we then have / = as 
expected). In the next section, we will show how the interaction range is cut off when the 
KAOSIS is below threshold. In the mean time though, we have already disposed of the reso- 
nant divergence associated with the process of Fig. ffl(b), which always has a KAOSIS above 
the one photon pair-production threshold. Indeed, we may now give a useful interpretation 
of that divergence: it arose because the intermediate photon (which at the KAOSIS may be 
viewed as being due to a pair annihilation) was not instructed to decay to a pair in a finite 
time, so that it produced finite amplitude for the final pair at all values of x 1 . The resulting 
total cross-section was infinite. Now that the dressed propagator is used to calculate the 
S-matrix element, the photon is aware of its decay obligations, and the cross-section due to 
this process is finite. 



V. CORRECTIONS TO SCATTERING STATES 

The propagator corrections described in the previous section are sufficient to control 
resonant divergences in many processes and regimes. Nevertheless, there are still cases where 
a scattering process may lead to a divergent resonance. Specifically, any process exhibiting 
a KAOSIS with vanishing decay rate will exhibit divergent scattering cross-sections even if 
calculated using the dressed propagators. 

Such processes are not at all rare. Consider the case of an electron scattering with 
a photon which is above the one photon pair-production threshold. The following two 
first-order processes correspond to an on-shell second-order process: first the photon pair- 
produces, then the initial electron annihilates with the newly produced positron to produce 
the final photon [see Fig. |2|(a)]. If the intermediate on-shell positron is in the Landau ground 
state, its decay rate is zero, so that there is no line width supplied by the propagator to 
control the resulting divergence. This divergence represents something of a calamity, since 
its effect is to make the total cross-section for Compton scattering divergent everywhere 
above the one photon pair-production threshold [|l7l ,|lS|1 . 

A second example is provided by electron-electron scattering, in which at least one of the 
initial electrons is in an excited Landau state [Fig. 0(b)]. Once again, there is a second-order 
on-shell process that is analogous to a succession of two first order processes, in which the 
excited electron emits a cyclotron photon which is in turn absorbed by the other electron. 
If the on-shell photon is below the one photon pair-production threshold, its decay rate is 



zero, and there is a resonant divergence — see the expressions in Appendix C of Ref. |22 . 

It turns out that in every second-order process with a stable KAOSIS there are non-zero 
decay rates associated with the initial and final scattering states. Thus, in the first of the 
two examples above, the initial and final photons are capable of decay, while in the second 
example, there are excited Landau levels in both the initial and final states. One might 
hope, then, that the decay rates of the scattering states might be pressed into service to 
control the resonant divergences when the decay rate of the intermediate state is zero. As 
we now demonstrate, this is not only possible, it is a necessary feature of the same program 
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of radiative corrections that brought the decay rates into the propagators. That program 
demands that we should apply radiative corrections to the external lines, in addition to the 
propagators. 

It might be objected at this point that loop corrections to external lines can have no 
bearing on the problem, since the arbitrary constants which arise during renormalization 
are in part fixed by the requirement that the external lines should suffer no corrections, so 
that the scattering states should continue to be represented by solutions of the "free" wave 
equation with the physical mass [[HJ. 

The answer to this objection is that in the present case, the limited number of arbitrary 
constants is not sufficient to satisfy the physical requirement above for all scattering states. 
For example, in the case of the electron field, we may only demand that the Landau ground 
state Ac propagate according to Eq. ( |2.1| ). Once we have used up the relevant renormaliza- 
tion constants to ensure this condition, the remaining excited Landau states must propagate 
according to Eq. ( ft.ip . In other words, in Eq. ( |3.1| ) we may set E(p° = r]E G , A Gl rf) = 
(where E G is the energy of the ground state), but then the remaining A ^ A G will yield 
non-zero on-shell, diagonal matrix elements of S. Similarly for the Maxwell field, we may 
only demand that II (u>k, k, j) = in the limit k — > 0, so that only infinitely soft photons see 
no refringence in the magnetized vacuum. 

The origin of this "feature" of magnetic QED is the fact that excited scattering states 
are technically not scattering states at all, insofar as they do not correspond to asymptotic 
one-particle states of the quantum field which are stable. In principle, we should only use 
stable states as initial and final scattering states — electrons and positrons in the Landau 
ground state and photons below the one photon pair-production threshold. The consequence 
of this restriction on the theory would be that multiple scattering events and scatterings 
followed by multiple emissions could only be treated by computing scattering amplitudes 
corresponding to very high order Feynman diagrams, a notoriously burdensome task. Thus, 
some of the most interesting astrophysical applications of the theory would become virtually 
inaccessible. 

As an alternative, we may represent such high-order multiple events as a succession of 
lower order transitions between states that are not necessarily stable, and treating those 
states as if they were genuine scattering states. For example, a process in which an electron 
in the ground state and a photon of energy above the third cyclotron harmonic make a 
resonant transition to a state with an electron in the ground state and four photons may be 
approximated by a Compton scattering event in which the final electron is left at the third 
harmonic, followed by three resonant decays. 

This approximation of scattering states by excited states is a common one in the literature 
P"8|^2T||, but to date there has been no investigation of its validity and limitations, or of what 
modifications the usual Feynman perturbation theory must suffer in order to accommodate 
them. That investigation is the central concern of this work. 

We now discuss the specific modifications to the scattering states due to radiative cor- 
rections. Consider a bare external electron line in a Feynman diagram which is represented 
in the scattering amplitude by the spinor a solution of Eq. (|2.1| ). After subjecting 

the line to corrections associated with the self-energy operator E, the result is a dressed line 
represented by the spinor Q^\x), where 
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= *5 ) + G fl .E.e^. (5.1) 

If we view the L-D operator [the wave operator in Eq. (|2.1]) 1 as a "free Hamiltonian" and E as 
a perturbation, Eq. (|5.1|) may be cast as a four-dimensional Lippmann-Schwinger equation 
connecting an eigenstate of the L-D operator with eigenvalue zero, to an eigenstate 0^ 
of the dressed L-D operator [the wave operator in Eq. (3.1)1, also with eigenvalue zero. In 
other words, 0^ satisfies 

[ 7 "(i0„ - eAj — m — E] ■ e<? = 0. (5.2) 

It is a simple matter to find solutions of this equation, since we already know of states that 
simultaneously diagonalize the "free" operator and its perturbation. Substituting 0^\x) = 
e- ip ° x ° '(j)^ (x) in Eq. we find 

p° =n{E A -iT A /2), (5.3) 

so that the dressed scattering state is 

ef } (i) = e-^-^^'VJ^x). (5.4) 

We may repeat the above argument for a bare external line which is represented by the 
Dirac conjugate spinor \l/^(x). The result is that the dressed state is A^\x), where 



Note that Ajp(x) is not the Dirac conjugate spinor of 0^ (x), since the real part of the 
exponent changes sign. 

The procedure is analogous for the external photon lines. The dressed states a(x) u are 
solutions of 

(g^D - IVK = 0. (5.6) 

Using the polarization states of the previous section (for j =_L, ||) to write 

a kj (xY = (2u k L")- 1 / 2 (e j re^ k ° x0+ ^, (5.7) 

and substituting in Eq. (|5.6j), we obtain 

k° = ±[u k -iT(k,j)/2], (5.8) 

so that the dressed scattering state is 

aiKj {x) v = (2cu fc L 3 )- 1 / 2 ( ej )"e ±i[ ^~ ir(kj ' )/2lx, °- ik - x . (5.9) 

From Eqs. ( |5.4j ), ( |5.5| ), and (|5.9|) it is apparent that the prescription of Eq. ( |1.1|) applies 
equally well to scattering states as to propagators. It is a general feature of this prescription 
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that the resulting positive and negative energy states are not conjugate to each other, since 
the real part of the exponent changes sign. Consequently, positive energy solutions decay 
as time increases, while negative energy solutions grow. This is in keeping with the inter- 
pretation of the negative energy solutions as particles which move backwards in time. In 
the application of these formulae to the calculation of S-matrix elements, the e ~ t ( E -' ir / 2 ) x 
dependence is ascribed to initial states, while the e + l ( E - lF / 2 ) x dependence is ascribed to 
final states. 

There is a second way of understanding the introduction of decay rates in the time- 
development exponentials of scattering states. We may take the view, discussed above, that 
the metastable scattering states are really approximations standing in for internal lines of 
much larger Feynman diagrams. That being the case, their functional form may be read 
directly from the components of their respective propagators, written in the forms of Eqs. 



(PD and (grp . 

The appearance of the decay rates in the time-development exponentials of the states 
will ultimately lead to their appearance in energy denominators of scattering amplitudes, 
after integration over time variables. Here, however, there appears a major difference with 
the usual practice of obtaining amplitudes. The real parts of the exponentials will lead in 
general to divergent expressions if the time integration limits are allowed to go to ±00 as 
usual. It is not difficult to see physically why we should expect trouble in this limit. Letting 
the upper time limit go to infinity in the S-matrix element is tantamount to asking the 
question, "what is the probability of observing this decaying state in the infinitely distant 
future?" Clearly, no calculation is required to see that the answer must be "zero." Similarly, 
letting the lower time limit go to minus infinity amounts to inquiring about an interaction 
at a finite time of a decaying state which was prepared in the infinitely distant past — a 
process which also has vanishing probability of occurrence. 

It is therefore necessary that the scattering theory be formulated for states prepared at 
finite times T; and measured at finite times Tf. Only when the initial state is stable may 
the limit Tj — > — 00 be taken, and only when the final, measured state is stable may we set 
Tf — > 00. These limits are termed the "absorption" and "emission" limits, for reasons that 
will shortly be made clear. 

Note that due to the real parts in the time-development exponentials, the scattering 
states are not generally normalized to unit probability. In fact, they may only be so nor- 
malized at a given, fixed time. Physically, it is necessary to ensure that the initial states 
be normalized at the preparation time T i: and that the final states be normalized at the 
measurement time Tf. This could be accomplished by setting x° —> x° — T^ or x° — > x° — Tf, 
as appropriate, in the expressions of Eqs. Q5.4|), (|5.5|), and (|5.9|). It is often simpler to cal- 
culate the amplitudes with the states as written above and multiply the result by the factor 
exp[|(TjTj — TfTf)], where ^i(f) is the sum of the decay rates of the particles in the initial 
(final) state. 

There is an extremely important consequence of this finite-time formulation of the the- 
ory: strict energy conservation no longer holds. The energy conserving 5-functions which 
appeared in the old amplitudes were a consequence of the integration of time exponentials 
over all time. By the time-energy uncertainty relation, we may not determine the energy 
to infinite precision over a finite time interval. This is not an alarming consequence of the 
theory, but rather a desirable one. When we discuss the excitation or decay of metastable 
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states, we cannot expect to determine the energy of those states to better than the natural 
line width, so energy-conserving 5-functions should actually violate our physical intuition for 
these processes. In fact, we will see that strict energy conservation is recovered only when 
Ti = Tf = 0. Thus, for example, in the expressions of Herold [|17]] for Compton scattering 



from ground state to ground state, the energy-conserving 5-functions are appropriate. 

We now discuss an example which provides a simple application of these ideas: cyclotron 
decay. Consider an electron which is prepared in an excited state A at a time Tj = 0. We 
calculate, to first order, the probability that the system should make a transition to the 
ground state Aq with the emission of a photon in the state (k, j), which is below the one 
photon pair-production threshold. Since Tf = 0, we may take the limit Tf — > oo. The 
S-matrix element is given by 



POO 

Sfi = T^ky e,-)" / dx° e -^ A ^T A /2)-E G -^ 
Jo 



-zT4 GA (ky ei r 

E A - E G - uj k - zTa/2' 
where T^ G ^(k) At is the interaction matrix element for the transition: 



(5.10) 



T AaA (kr = e 2 L" 3 / 2 (2^)- 1 / 2 J ,f 3 x e - ik ^W(x) T ^ +) (x). (5.11) 

Note that in Eq. ( |5.10| ), the energy conserving 5-function has been replaced by the Wigner- 
Weisskopf line shape function. Thus, the new formalism has reproduced the well-known 
result from non-relativistic quantum mechanics, which the old formalism could not (compare 
Refs. I-Pl!). 

We close this section with a discussion of the effect of using the dressed electron scattering 
states of Eqs. (|5.4j) and (|5.5|) on the range of the interaction in the presence of a KAOSIS of 
the photon propagator below the one photon pair-production threshold. We may attempt to 
repeat the procedure that led to Eqs. ( |2.5| ) and fl4.13Q . However, this time the computation 
of the S-matrix element corresponding to Fig. (Tj(a) is no longer equivalent to taking the 
Fourier transform of the photon propagator with respect to x°, since the dependence of the 
scattering states on x° is no longer purely oscillatory. Rather, the effective one-dimensional 
interaction g(x l ; k 2 , k 3 ) that results is given by 

/dk 
— hix'-^^k^k^W^-E), (5.12) 

where W(z) is a function which is only appreciable in a range ±A about z = 0. The 
restricted domain of W is a reflection of the restricted domain in time of the scattering 
states — in fact, we have either A ~ T or A ~ (Tf — Tj) _1 , whichever is largest. 

By using stationary phase arguments, it is easy to see that g(x 1 \ k 2 , k 3 ) can only be 
appreciable for a limited range of | rc 1 1 : 



2tt 



\x 



i| < (5 13) 

~ [(E + A) 2 - (P) 2 - (P) 2 ] 1/2 - [£ 2 - (A; 2 ) 2 - (P) 2 ] 1/2 ' 
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Thus the range of the interaction is curtailed when the scattering states may decay. We see 
that the spurious infinite interaction range that entered calculations which used bare excited 
scattering states was a consequence of their assumed infinite duration. Once their limited 
duration is incorporated into the formalism, their interaction becomes short-ranged. 

Note that a KAOSIS of the photon propagator may exist only if either some of the 
scattering states are excited or the KAOSIS itself is above the one photon pair-production 
threshold. Therefore, the scattering cross-sections for processes with virtual photons are 
now always finite if the dressed states and propagators are used to calculate them. We will 
show in the next section that resonant divergences are now also under control in processes 
with virtual fermions. 



VI. SCATTERING FORMULAE 

We now compute generic second-order formulae for two-particle to two-particle scatter- 
ing. The Feynman diagram in Fig. depicts such a process irrespective of whether the 
various lines represent fermions or photons. Our notation, which is illustrated in Fig. |3|, 
is as follows. Let the energies of the external lines be E p , and let their decay rates be T p 
(p = a, b, c, d). Define p p = ±1, with p p = +1 if the line represents an incoming state and 
fi p = —1 if it represents an outgoing state. Let the lines with p = a,b join at the vertex 
with coordinates x, and those with p = c,d join at the vertex with coordinates y (see Fig. 
|3]). Define E x = p a E a + PbEb, T x = /i a r a + pb^b, and similarly for E y and Y y . Also define 
the complex energies £ x = E x — ^T x , £ y = E y — ^T y . 

Let the energies of the two initial particles be en, e^, and let their decay rates be jn, 
7i2- Also let the energies of the two final particles be e/i, e/2, and let their decay rates be 
7/i> If 2- The e and 7 are set equal to the E and Y as appropriate to the process under 
consideration, an identification illustrated by the passage from Fig. |3] to Fig. f|. 

We denote the quantum state of the intermediate particle by I: I = A if the particle is a 
fermion, I = (k, j) if it is a photon. The energy and decay rate of the intermediate state are 
Ei and Ti, respectively, and we define £1 = Ei — We also introduce the index r] = ±1, 
where r] = +1 if the intermediate state has positive energy and rj = — 1 otherwise. 

The propagators are given by Eqs. ( p.5|) and (|4.12| ). The S-matrix element for the process 



is obtained by sandwiching the appropriate propagator between the scattering states in the 
usual way [Q and integrating over the space time coordinates x and y. The general result 



is 

S fl = tJ2c v M fl (l, v )^$ f (£ x ,£ y ). (6.1) 

Ir, 

Here, Mfi(l, rj) is the product of two interaction matrix elements appropriate to the process, 
c v is 1 if the intermediate state is a photon and —7] if it is a fermion, and 



AlV) (£ g ) = e +hii+~fi2)Ti/2-{-y fl +j f2 )T f /2 +i(e il +e i2 )T i -i(e fl +e f2 )T f 

dx° [ 1 dy°9 [rj (x° - y )] exp [-vqE^x - y°) - i£ x x° - i£ y y°] . (6.2) 



x 



Ti JT, 
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The factor e + ( 7il+7i2 ) Ti / 2 (t/i+w) t // 2 in Eq. (|6.2| ) adjusts the normalization of the initial and 
final states to be 1 at Tj and T/, respectively, while the factor e + ^ en+ei2 ' )Tl ~' l ^ efl+e f^ Tf allows 
a convenient choice of phase. The quantity ^j,.{£ x , £ y ) is the new object which incorporates 
the resonant energy denominators and in general replaces the energy-conserving 5-function. 
It may be calculated by the substitution u = x° — y°, v = x° + y°. The result is 

(n +(7il+7i2)T i /2-(7 / i+7 /2 )T / /2„+i(e i i+ei2)T i -i(e / i+e /2 )r / 

£v)= 

e -i(£ x +£ y )T f j e (i/2)[(l- J7 )£: !c +(l+T,)£: y -2£ ! ](T / -T i ) _ jl 

(1/2) [(1- rj)£ x + (1 + rj)£ y - 2£ t ] 

e -i(£ x +£ y )Ti ^ e (i/2)l(~l- V )£ x + (-l+v)£y'^im-T t ) _ ]] 

(l/2)[(-l-ri)S x + (-l + ri)S y -2S l \ 



(6.3) 

Eq. (|6.3|) provides an expression which may be adapted to any second-order process of 
interest by adjusting the /i and the assignments of the e and 7 to the E and T. In the 
particular case of two-particle to two-particle scattering we need consider four cases: the 
"direct" and "cross" diagrams (Fig. |j), each for 77 = ±1. The expressions are simplified by 
the notation Ae = e a + e i2 - e fx - e f2 , A7 = j a + j i2 - 7/1 - 7/2, and r = T f - T t . 

a. Direct diagram, 77 = +1 We obtain this case from Eq. ( |6.3| ) by setting £ x = — (e/i + 
e/2) + 2(7/1 + 7/2) /2 and £ y = (e n + e i2 ) - 2(7*1 + 7*2) /2. We obtain 

f {,+ ) (£ £ ) = x <^ - - 

«T<2> \ x: ^ Ae _ zA7/2 I (eii + ^ _ Sj) _ i{rfii + ^ _ rj)/2 
p -i{Ei-iTi/2)T _ p -t[(e/i+e/a)-i(7/i+7/2)/2]7 



(e/i + e/ 2 - £2) - 2(7/1 + 7/2 - 1/0/2 



(6.4) 



6. Direct diagram, rj — — 1 For the same assignments as case (§), but with 77 = —1, we 
find 



'VTiTf {£x,£y) — Ag _ 



iA 7 /2 

e -*[(ea+ei2+e/i+e/2+-E ; ;)-*(7ii+7i2+7/i+7/2+ri)/2]T _ g— i[(eii+e<a)— »(7il+7i2)/2]T 

(-e 7 i - e/2 - -Si) - 2'( _ 7/i - 7/2 - r z )/2 

e -i[(eii+ei 2 +e/i+e / 2+S i )-i(7a+7i2+7/i+7/2+r ; )/2] _ g-i[(e/i+e/2)-*(7/l+7/a)/2]7 

(-eji - e i2 - Ei) - i(-7ii - 7 i2 - r t )/2 



(6.5) 



c. Cross diagram, 77 = +1 Here we set ^ = (e^ — e/2) — i(7«2 — 7/2)/2 and £ y 
(eii - e/i) - z(7a - 7/i)/2: 



,, I f g-i[(e/i+ei 2 +-E0- i (7/l+7i2+r ! )/2]r _ „-i[(e il +ei2)-i(7ii+7i2)/2]T 

^T f V*, ty) Ae _ zA7/2 { ( g a _ e/i _ ^) _ ..^ _ _ rj)/2 

-*[(e/i+ei 2 +B()- i (7/i+7i2+ri)/2]r _ R -i[(e /1 +e /2 )-«(7/l+T/2)]7" 

(6.6) 



(e/2 - e i2 - £/) - 2(7/2 - 7i2 - r z )/2 
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d. Cross diagram, 77 = — 1 The assignments here are as in case and 77 = —1: 

,, 1 f p -«[(e/2+eii+-B i )-j(7/2+7ii+r i )/2]T _ -i[(e l i+e i2 )-«(7 I i+7i2)/2]7 

rCi-iJ/c- <? \ _ _ w J _ 



$w( e *> e v)- Ae-iA T /2 X { 



(e<2 - e/ 2 - Ei) - i(<y i2 - 7/2 - r,)/2 

e -«[(e/2+eii+-B;)-*(7/2+7 l i+rO/2]T _ e -i[(e/i+e /2 +,E i )-i(7/l+7/2+r ! )/2]T 

(e/i - e a - - £(7/1 - 7a - r { )/2 



(6.7) 

The expressions in Eqs. ( |6.4|) , ( |6.5|) , ( |Q>D , and (|6.7|) all contain products of two complex 
energy denominators, one "energy-conserving" and the other resonant. All these energy 
denominators may in principle go to zero. When this happens, however, there is always 
a cancellation of the exponentials in the numerators, so that the result is always finite. 
This is a consequence of the fact that these expressions were derived starting from Eq. 
( p.2|) , an integral over a finite range of a finite integrand, which consequently may never 
diverge. Therefore, scattering processes containing virtual fermions now have finite S-matrix 
elements. These were the last remaining resonant divergences in the theory, which is now 
entirely finite. 

We have thus eliminated all resonant divergences from our scattering cross-sections. The 
price we have paid is the dependence of the S-matrix elements on the time lapse r between 
the preparation of the initial state and the measurement of the final state, and the attendant 
loss of strict energy conservation. Note that in general, the nature of the dependence on r 
is for the S-matrix elements to decay away as r — > 00. As discussed in the previous section, 
this is the behavior expected for scattering from excited states to excited states. 

Note also that these expressions lack crossing symmetry. The reason is the introduction 
of the exponential factor outside the integral in Eq. ( |6.3|) , which is not symmetric under 
the replacement e»i «-> —e/2, ju <-> —7/2- If we divide the expressions of Eqs. (|6.4|), (|6.5|), 
( p.6|) , and (|6.7| ) by the exponential factor we find that the resulting expressions are, in fact, 
crossing symmetric. 

In order to parlay the above expressions into cross-sections, we must substitute them into 
Eq. ( |6.1D to obtain a r-dependent expression for the S-matrix element. The reaction rate R 
is then given by R = d\Sfi\ 2 /dr = 2Re{Sfi*dSfi/dr}. The cross-section may be obtained 
from R by the usual kinematic manipulation: da/dQf = L 3 |v| _1 i2, where dQf is a volume 
element in the space of final states and v is the relative velocity of the initial particles. 

There are two special cases where it is possible to eliminate the dependence on r from 
these expressions: when the initial state is stable (7^ = 7^ = 0), and when the final state 
is stable (7/1 = 7/2 = 0). These cases are called "absorption" and "emission" scattering, 
respectively. In the absorption scattering case, we may set T{ — > —00, while in the emission 
scattering case we may set 7} — > 00. In either case, r — ► 00, and the above expressions for 
the £ are greatly simplified. 



A. Absorption scattering 

The special case of a stable initial state corresponds to a situation analogous to absorp- 
tion, in which we prepare beams of stable particles and observe the excited products before 
they have the opportunity to decay. If we set ja = 7i2 = in Eqs. ( |6.4|) , (|6.5|), ( |6.6|) , and 
(p.7|), and take the limit r — > 00, we find the following results: 
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For the direct diagram, 



■-"^ Ey) ~ (Ae - iAj/2)[(-e n - e l2 + E t ) - ' (6 ' 8a) 



tt>%&> Q (Ae _ iA 7 /2)[( e/1 + e/2 + e,) - i( 7/1 + 7/2 + r,)/2] ' (6 ' 8b) 

and for the cross diagram, 



d'S^Sy) - {Ae _ iAl/2)[{efl _ ea+El) _ i{lfl + ri)/2 y (6 - 9a) 



£y) - (Ae _ ,A 7 /2)[(e /2 - e i2 + E l ) - i{ lf2 + r,)/2] ' (6 ' 9b) 

We have eliminated an inessential phase factor. Note that A 7 = — ( 7 /i + 7/2)- Substituting 
these expressions in Eq. fl6.1| ), we obtain a time-independent S-matrix element. In order 
to get the result into the form of a cross-section, note that the final state is decaying at a 
rate 7/1 + 7/2- In order that the probability of that final state be time-independent and 
equal to |5/i| 2 , the reaction rate must exactly balance the decay rate, so that we must have 
R = (7/1 + 7/2)|S/i| 2 = |A7||S/j| 2 . Again, cross-section da/dQf may be trivially obtained 
from R. Note that da/dQf contains the following functional dependence on Ae: 

da IA7I A7-^o 

oc — j--r — ► 27rd{e il + e i2 - e fl - e /2 ). (6.10) 



dn f Ae 2 + (A7/2) 2 

In other words, the energy conservation in the cross-section is Lorentzian, and in the limit 
of stable final states we recover the energy-conserving 5-function with the correct coefficient 
of 27T. As expected, the strict energy conservation in the old expressions for the S-matrix 
elements is correct for scattering from stable states to stable states. 



B. Emission scattering 

The special case of a stable final state corresponds to a situation analogous to emission, 
in which we prepare a beam of excited particles and observe them after they have scattered 
into the stable final state. We set 7/1 = 7/2 = in Eqs. (|6.4|), ( |6.5| ), (|6T^|), and (|6.7|), and 
take the limit r — > 00, to obtain: 

For the direct diagram, 



egj?^ Sy) - (Ag _ ,A 7/ 2)[( e/1 + e f2 + E t ) - ^/2] ' (6 ' lla) 
&J{£*,£y) = (Ae _ ,A 7 /2)[(e 4l + e,, 2 + E t ) - ifa + 7,2 + r,)/2] ' (6 ' llb) 
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and for the cross diagram, 

Z%£\£;£v) = ( Ae - *A 7 /2)[(e i2 - e /2 + E t ) - *( 7i2 + r,)/2] ' ( ° 2a) 



^oo 1 - (Ae _ iAl /2)[(e a - e fl + E l ) - i( 7tt + r,)/2]' (6 ' 12b) 

Once again, we have eliminated an inessential phase from the amplitudes. Now we have A7 = 
la + 7i2- Substituting these expressions in Eq. Q6.1|) , we again obtain a time-independent 
S-matrix element. We obtain a cross-section by the argument illustrated in Fig. [5], which 
depicts a semi-infinite tube of cross-sectional area da/dQf, terminating at the position of 
the target particle and extending in the direction of —v. The total probability (per unit 
final phase-space volume) of an interaction leading to a final state in dQf is equal to the 
integral over the interior of the tube of the probability that each infinitesimal slice should 
contain the projectile particle and that it should actually reach the target particle in spite 
of the fact that the two-particle state is decaying away: 

\b J% \ — L J |v|<tt e -A 7 L3diy {bA6} 



so that 



<h L 3 \v\- 1 A 1 \S fi \ 2 . (6.14) 



dn 



We again see the Lorentzian energy conservation of Eq. ( |6.10D , so that in this case we 



also recover strict energy conservation for stable state to stable state scattering. 



VII. DISCUSSION 

The regime of validity of the "emission" and "absorption" scattering limits is obvious 
from their context. On the other hand, the general formula in Eq. ( |6.3| ) requires some 
discussion. As discussed previously, the formula never "misbehaves", in the sense that it 
never yields a divergent result. In fact, in general that result tends to zero as r — > 00. 
While this is physically sensible, it obviously makes the large time-lapse limit a less than 
useful one. It is clear that what fails in this limit is the validity of the perturbation-theoretic 
order of the calculation. Since the scattering states are themselves decaying to other states, 
those other states should be included in the calculation, leading to higher-order processes. 
The second-order calculations outlined in the previous section are only useful for values of 
r such that the scattering states have little chance to decay, that is for IV <C 1, where T is 
the largest of the decay rates in the process. For example, we might choose r to be on the 
order of a collision time, if we are studying a gas with density and temperature such that 
the collision rates far exceed the decay rates. 

The general case above obviously represents a fairly radical departure from the usual scat- 
tering formulae. The "emission" and "absorption" scattering limits amount to somewhat less 
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radical modifications. One obvious such modification is the replacement of strict energy con- 
servation with "Lorentzian" energy conservation. Another is that the "non-resonant" energy 
denominators [Eqs. ( |6.8b| ), ( |6.9aj ), ( |6.9b| ), ( |6.11b| ), ( |6.12a| ), ( |6.12b| )1 now contain the decay 
rates of the scattering states as well as those of the intermediate states. The importance 
of these changes depends upon whether the decay widths entering the energy denominators 
are electron decay widths or photon decay widths. 

If the decay widths are purely fermionic, they are gently varying functions of energy, 
and their magnitudes are smaller by e 2 than their own characteristic scale of variation, the 
scale of variation of the interaction matrix elements, and the characteristic separation of 
the resonances. Consequently, in this case the relative change that results from introducing 
Lorentzian, rather than exact, energy conservation, and from introducing the decay widths 
of the scattering states into the energy denominators, is of order e 2 . On the other hand, 
if some of the decay widths correspond to photon lines, they can vary rather rapidly as 
a function of energy [|TT|,|T^]. Thus, the behavior of the energy denominators is not really 
"Lorentzian" , despite notational appearances to the contrary. The departure from the cross- 
sections computed assuming strict energy conservation and not including external line decay 
widths might turn out to be appreciable in this case, although its precise magnitude remains 
to be assessed. 

In the limit of stable scattering states the usual results are completely reproduced, since 
we recover strict energy conservation and there are no scattering state decay widths to 
include in resonant energy denominators. 

In connection with the dressed photon propagator of Eq. (|4.11|) , we wish to comment on 
a point which is a potential source of confusion. The decay width r(k, j) is to be evaluated 
on the light cone, as implied by the first line of Eq. ( 4.10 ). Now, when the photon propagator 
is used in an S-matrix element, the values of k°, k 2 , and k 3 are fixed by the x°, x 2 , and 
x 3 momenta of the scattering states. Thus, the sum over intermediate photon wave states 
involves an integral over the component k l of the wave vector k. The value of T(k, j) must 
be evaluated for each value of k 1 in the integral. This is analogous to the case of the electron 
propagator, Eq. ( |3.4|) , in which the intermediate electron decay width is evaluated, on the 
energy shell, for each Landau level in the sum over intermediate states. The only difference 
between the two cases is that the relevant degrees of freedom are discrete for the electron 
propagator, while they are continuous for the photon propagator. 

The radiatively corrected photon propagator permits for the first time the evaluation of 
processes such as e + e~ — > e + e~, e~ — > e~e + e~, and 7e~ — > e~e + e~, all of which are important 
for neutron star emission. Of equal astrophysical importance is the evaluation above the 
one photon pair-production threshold of Compton scattering, two photon pair annihilation, 
and two photon pair production, which is now possible by virtue of the radiatively corrected 
scattering states. Finally, we now have access to the processes e~e~ — > e~e~ and e + e _ — > 
e + e~ even when the initial and final states are excited. 
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FIGURES 



FIG. 1. Two diagrams for e + -e~ scattering. The labels a+, denote the x -coordinate of the 
orbit center of the initial positron and electron, respectively, while the labels , a J denote the 
orbit center of the final positron and electron, respectively. The spatial separation between the two 
fermion lines is given by the parameter s = — a J for (a), and by t = a J — for (b). In either 
case the interaction does not fall off with increasing separation, so that the result is a divergence 
in the total cross-section. 

FIG. 2. Processes allowing non-decaying, on-shell virtual states; (a) e~7 scattering; (b) e~e~ 
scattering. The process in (a) may be viewed as a pair production followed by a pair annihilation. 
When the intermediate positron is in the Landau ground state it has zero decay width. Similarly, 
the process in (b) may be viewed as cyclotron emission followed by cyclotron absorption. The decay 
width of the intermediate state vanishes if it is below the one photon pair-production threshold. 

FIG. 3. Generic scattering diagram. The lines may be either fermion or photon lines, and the 
external lines may be either incoming (fi = +1) or outgoing (fj, = — 1). 

FIG. 4. Generic "direct" and "cross" diagrams. The lines may be either fermion or photon lines. 

FIG. 5. Kinematics of emission scattering. The target particle is at the end face of the 
semi-infinite tube of cross-sectional area da/dUf. The projectile particle must be in the tube 
in order for an interaction to occur. The decreasing density of circular sections in the figure is 
meant to represent the decreasing probability of an interaction due to the exponential decay of the 
initial state. 
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